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Abstract 


It was recently shown that there exists nretastable U(1 )a domain wall con¬ 
figurations in high-density QCD (/j > 1 GeV). We will assess the stability 
of such nontrivial field configurations at intermediate densities (p < 1 GeV). 
The existence of such configurations at intermediate densities could have in¬ 
teresting consequences for the physics of neutron stars where such densities 
are realized. 


I. INTRODUCTION 

In general, there do not exist domain walls or other topological defects within the Stan¬ 
dard Model. This is due to the trivial topology of the vacuum manifold. In contrast, these 
objects are quite common in condensed matter physics and cosmology. However, it has been 
realized only recently jlj-Q that topological defects such as domain walls and strings may 
exist within the Standard Model at a large chemical potential. 

It is well known that topological defects result from symmetries being broken. In the past 
few years there has been renewed interest in high density QCD. Similar to the BCS pairing 
in conventional superconductivity, the ground state of QCD at high density is unstable due 
to the formation of a diquark condensate 0-0 (see H for a review). In this new ground 
state, various symmetries which are present at p = 0 are broken by the presence of this 
non-zero diquark condensate. This leads to the formation of the various topological defects 
discussed in [|1]-Q]. 

In [|TJ, it was shown that at high densities (p 1 GeV) there exist domain wall solutions 
which interpolate between the same vacuum state, in which the U(1 )a phase of the diquark 
condensate varies between 0 and 2n. This type of domain wall which interpolates between 
the same vacuum state has been studied before in the context of axion models ||. It is 
interesting to note that similar domain wall configurations are present for the zero density 
case in the large N c limit and could be present for the physically relevant case of N c = 3 |IU| , 
in which case they can be studied at the BNL Relativistic Heavy Ion Collider (RHIC) ||11|| . 
Given this, one might ask what happens between these two regions of p = 0 and /i> 1 GeV. 
The main goal of the present paper is the analysis of the classical stability of the U(1 )a 
domain walls for p < 1 GeV. We demonstrate that the U(1 )a domain walls are classically 
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stable down to densities /i ~ 800 MeV. The main idea is to interpolate between 1 GeV 
(where perturbation theory is justified) and /i ~ /i c (where instanton calculations lead to a 
reasonable description). Ideally, we would like to be able to make definitive statements on 
the stability of such configurations all the way down from large /i to /i c ~ 500 MeV, where 
/i c is the critical chemical potential above which the color superconducting phase occurs. 
Unfortunately, we are unable to make a definitive statement for /i c ^ n < 800 MeV due to 
a lack of theoretical control in this region. One can only speculate on the behavior in this 
region. 

This work is organized as follows. In Sec. 2, we will construct the effective potential. In 
Sec. 3, we will describe the nontrivial domain wall solutions. In Sec. 4, we will examine the 
classical stability of such configurations under small perturbations. In Sec. 5, we will end 
with concluding remarks and future considerations. 


II. THE EFFECTIVE POTENTIAL 

It is well known that the ground state of Nf = 2,3, N c = 3 QCD exhibits the Cooper 
pairing phenomenon as in conventional superconductivity (known as the color 2SC (2 flavor 
color superconducting) and CFL (color-flavor locked) phases of QCD B0.0H)- In what 
follows, we will consider the N c — Nf — 3 CFL phase. In the CFL phase, the condensates 
take the form 


, 4 >* ~ 

(&«&>' ~ (i) 

where L and R represent left and right handed quarks; a, /?, and 7 are the flavor indices; i 
and j are spinor indices; a, b, and c are color indices; and A " 7 and V c 7 are the condensates 
which are complex color-flavor matrices. As the magnitudes of these condensates depend on 
the color index c, one can easily see that these objects are not gauge invariant by themselves. 
In order to construct a gauge invariant held, the following matrix which describes the octet 
of mesons and axial singlet was considered in [|l^,|T^]: 

j$ = xy' = '£xI!y?. ( 2 ) 

c 

If a U(1 )a rotation (q —> exp * 75 °/ 2 q) of this gauge-invariant held £ is performed, we see that 
the helds (|l|) transform as 


X e~ ia X, 

Y -»■ e +ia Y, (3) 

and therefore 

£ e~ 2ia Y. (4) 

Goldstone’s theorem states that there must be a single Goldstone mode rf associated with 
the breaking of this symmetry. Given this, we can parametrize the held as follows: 
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E = E 0 e p e~ l(p , 


( 5 ) 


where the phase 0 — v'/ f is defined as a dimensionless held which describes the rf boson, / is 
the corresponding decay constant, E 0 is the vacuum expectation value of the composite held 
(||), and p is another dimensionless held which describes the fluctuations of the magnitude 
of the condensate^ (analogous to the a held related to the huctuations of the ( qq ) chiral 
condensate). We choose to parametrize the held as e p for convenience later on. 

In order to construct an effective potential describing the dynamics of the phase of 
the condensate 0 as well as the magnitude |E|, there are two types of terms which must be 
included. The hrst term which explicitly breaks the U( 1) J 4 symmetry was calculated in |l, 
by substituting the form of condensates given above into the instanton induced four-fermion 
Lagrangian [|T7| , |T^1 : 
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U-tort (/>,<!>) = e p cos 0 , 


( 6 ) 


where A 0 is the value of the gap in the quark spectrum. The perturbative form of the 
expectation value of the condensate has been used in arriving at this result [SLIT - 


|£|} = £« = 


9 


8n 2 g 2 

For Nf = 3, the dimensionless coefficient a was found to be 


( 7 ) 


a — 7 x i° 3 DA (ln^ 




QCD 




( 8 ) 


where m s is the mass of the strange quark. The mass of the corresponding rf boson can 
be easily calculated by expanding the potential, my = 27r-y/a A 0 . According to [|lj, this 
potential is only under theoretical control when the mass of the rf boson is much less than 
the typical scale for higher excitations 2A, which corresponds to a <C 1 / 7 T 2 . For physical 
values of the strange quark, this corresponds to a chemical potential of about /i ~ 700 MeV 
(for a ~ 1 / 7 T 2 ). 

To be able to draw any conclusions about the stability of the domain wall solution 
for intermediate densities (/i c ^ /i ^ 1 GeV), one must include degrees of freedom which 
are related to the fluctuations of the absolute value of |E|. We do not know the effective 
potential in the region of interest; however, for qualitative discussions we shall use a potential 
derived for asymptotically large p. The second type of term which must be included in an 
effective potential description is one which uniquely fixes the magnitude of the condensate. 
An effective potential for the magnitude of the condensate |E| was derived in 


in the 

perturbative region where the analytical form of the gap is known. This effective potential 
fixes uniquely the value of the vacuum condensate. The potential is of the Coleman-Weinberg 
type B and is given as follows: 


1 The field p should not be confused with the familiar p meson in QCD. 
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fi 
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7T 2 
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9tt ’ 


S 0 


e'U - p), 


(9) 


where a s is the standard strong coupling constant and the perturbative result of the con¬ 
densate (|7|) has been used. In all calculations that follow we will assume A 0 = 100 MeV as 
the numerical value for the gap. It should be noted that the potential ([]) is justified only 
in the region 


i'log(|E|/E„) « 1. 


( 10 ) 


Since we are considering E « E 0 , the use of this effective potential is justified. 

We are interested in the region p c < p < 1 GeV, where Eq. (|9|) is not literally correct. 
Even though the results stated above are not necessarily under theoretical control for p > p c , 
one can speculate on how the coefficients behave at intermediate densities when p c < p < 
1 GeV. Due to the fact that all the same symmetries are present as the chemical potential 
is lowered until reaching p c , we expect the qualitative form of the effective potential (,5j) 
and (^|) to remain the same. As the chemical potential is lowered, eventually perturbative 
calculations which are valid at asymptotically large p are no longer the correct description 
and instanton calculations become relevant. One would expect that these calculations must 
match up at some point. However, the coefficients in front of the potential could possibly 
be very different in the density region of interest. We will refer to the coefficient in front of 
the one-instanton potential as f3\ and the coefficient in front of the perturbative potential as 
0 2 . The value of is essentially fixed by the form of the condensate and by the constituent 
quark mass. Below the critical chemical potential p c at which the chiral phase transition 
occurs, this coefficient is independent of p |H§. Therefore, we would expect that 0i would be 
some smooth function of // which reaches its maximum value at /i — p c . ht p c < p < 1 GeV, 
the coefficient /3 2 is modified by the formation of instanton-antiinstanton (II) molecules |I3| . 
We should also note that the coefficient 0! can also be estimated from the instanton liquid 
model JHJ using the average size of instantons as well as requiring a constituent quark mass 
of about 350 — 400 MeV. 

Combining both terms we have an effective potential which is given by the following: 


V(p, 0) = -/3 2 e p ( 1 - p) - (3\e p cos 0. 


( 11 ) 


The values of the coefficients (3\ and 02 are known for asymptotically large p: 

( 12 ) 

(13) 

In Fig. |Tj, we show the cross section of the potential at 0 = 0, 27t and 0 = 7r for p = 800 MeV. 
Notice that the existence of an absolute minima at 0 = 0, 2ir and a saddle point at 0 = 7r 
allows for nontrivial configurations which wind around the barrier at |E| = 0. In the limit 


0i = ap 2 Al 


02 = 


7T 
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p —» oo, the parameter a —> 0 and the potential has degenerate minima at E = E 0 . The 
kinetic term is given by 


|<9o£| 

2 , 


U 2 <9;£ 2 


£ 

2 


(d oP ) 2 - u 2 (d iP ) 2 + (<9 O 0) 2 - « 2 (<9,0) 2 , 


(14) 


where u is the velocity which is different from 1. 
constant / and velocity u were calculated in 
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The perturbative values for the decay 
In order to fix the correct dimensionality, 
we must multiply the kinetic term by the appropriate powers of the decay constant. The 
full effective Lagrangian up to two derivatives in the fields is then given by 


C = f 2 [{d 0 p) 2 - u 2 (dip ) 2 ] + / 2 [(<9 O 0) 2 - u 2 (di(f)) 2 } - V(p, 0). (15) 


In the above we have assumed that fp — f p — f. The exact numerical value for f p is not 
known. However, in the large p limit fp ~ f p ~ p in order to have an appropriate scale for 
77i p ~ A 0 (once again, the p field should not be confused with the well known p meson in 
QCD at p — 0). 


III. DOMAIN WALL SOLUTIONS 

As was done in [|l|, if we replace the field |E| = £ 0 e p by its vacuum expectation value 
(which is justified for p 1 GeV), the resulting potential is of the sine-Gordon type. The 
Lagrangian is given by the following: 

C = / 2 [(<9 O 0) 2 - U 2 (di(t)) 2 } - Vl-inst (0), (16) 

where the constant term has been dropped. The static domain wall solution to the corre¬ 
sponding equation of motion is well known. Considering a domain wall in the z direction, 
the solution is given by 


4>o(z) = 4 arctan (exp(— mz/u )), (17) 

where m is the mass of the p'. This solution interpolates between the same vacuum state; 
at z = Too we have 0 = 0, 27r, respectively. It is well known that this solution is absolutely 
stable under small perturbations 0 = 0 o + 50. In other words, the Schrodinger-type equation 
obtained by varying the field and linearizing the equation of motion, 

- d 2 Scfrn + — (^ 1-2 sech 2 ^ zjj <50 n = u; 2 50 n , (18) 

has the lowest eigenvalue u> 0 = 0 corresponding to 5(p a = d<fi 0 (z)/dz ~ sech(m^/n). This is 
just the zero mode which is a result of translational invariance z —» z + z D . Since the lowest 
eigenvalue is non-negative, the domain wall solution is stable under small perturbations. It 
turns out that this is the only bound state which is a solution of Eq. ( |I8|) . 

In the case in which the replacement |E| —> E 0 is not done, the solution must be modified. 
If we want to study a stable solution for p which is not asymptotically large, we must include 
fluctuations in 0 as well as in the p direction (i.e., the absolute value of E). The two equations 
of motion for static solutions are given by 
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2 pu 2 V 2 p = p e p — Pi e p cosp, 

2 pu 2 V 2 p = Pi e p sin p. 


(19) 

( 20 ) 


Although we do not know the exact solution for this set of coupled nonlinear differential 
equations, if P 1 /P 2 < 1 we can approximate the solutions. In this case, the approximate 
solutions can be parametrized by 

e po m 1 + a cosp Q , (21) 

p a ~ 4 arctan (exp(— mz/u )), (22) 

where a ~ Pi/Pz- Our stability analysis will be based upon these approximate solutions of 
the equations of motion. Since the above solutions do not correspond to the exact solutions 
to the equations of motion (minimum energy path which winds around that barrier at 
|E| = 0), there will be nonzero linear terms when the energy of the system is perturbed 
about the domain wall solutions p D and p 0 . These will be estimated in the following section 
where the stability analysis is performed. 

IV. STABILITY ANALYSIS OF DOMAIN WALLS 

Although these domain walls may exist as classically stable objects at large densities, it 
is not immediately obvious if this is the case at intermediate densities. In order to examine 
the classical stability of the domain wall configurations in the region p < 1 GeV, we must 
look at how the system reacts to small perturbations of the fields. The energy density is 
given by the following expression: 

/ +oo 

dz[fu 2 {Vp) 2 + /V(V0) 2 + V{<j>,p)\. (23) 

-oo 

In order to express all integrals as dimensionless quantities, we will perform the following 
change of variables: 


2 = 


A’ 


A = 


u 


m v > 


The energy density is now given by 

f 2 u 2 


E(<j>, p) — 


A 


J &'((Vp) 2 + (V0) 2 + Pv(4,,p) 


f- 


(24) 


(25) 


where the derivative is now taken with respect to the dimensionless coordinate z'. We can 
see that this has the correct dimensions of MeV 3 times the dimensionless integral in square 
brackets. Following the standard method for analyzing the stability of a classical solution 
which was briefly described in the previous section, we will expand the fields about their 
vacuum expectation values: 


P -»• Po + dp, 

P —> po + dp. (26) 
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Next, we substitute Eq. (^ 6 [) into Eq. (|23| ) and perform an expansion, neglecting any terms 
greater that quadratic order in 8 p and 8 (f)-. 


E ~ £ (0) + E {1) + E (2) 


— E((f) 0 ,p 0 ) +7 J dz' 8 p\ 

+ 7 J dz' 8 (f) 

-2 V 2 (f ) 0 + 

+ 7 J dz 8 <f) 

( A 2 8 2 V 

1 f 2 u 2 8 (f> 8 p 


A 2 5E 


f 2 u 2 8 p 


<po,po / 


A 2 8 V 


<j>o,po / 


f 2 u 2 8 (f) 

8p + 7 J dz' 8 (f> 


+ 7 J dz' 8p 


-V 2 + 


A 2 8 2 V 


2 f 2 u 2 8 p 2 
A 2 8 2 V 


(j>o iPo / 


4>o,po 


2 f 2 u 2 8(f) 2 


8 p 

(27) 


<fioipo / 


where 7 = f 2 u 2 / A. The first term E ^ in the above expansion is the energy density or 
wall tension of the domain wall. In the case in which the domain wall solutions given by 
Eq. (|2T1) were the exact solutions to the classical equations of motion, the linear terms 
(proportional to 8 p and 8 (f> ) would be zero everywhere. Due to the fact that our solutions 
are not exact, these must be considered. 

First, let us estimate the term which is linear in 8 (f>: 

f 2 u 2 r ( A 2 A 

—J dz'8(f> (-2V 2 0 O + -j^(8ie po sin 0 o j . (28) 



We know that e Po ~ 1 + p Q and using the fact that <f ) 0 is a solution to the equation of motion 
given by Eq. (j2(]) with p = 0, we have: 

8 (f) {1] « f dz' /?ip o sin0 o 50 

A J j 2 u 2 

f 2 u 2 r 

~ J dz' a 7 r 2 sin 2 cf > 0 8 (f). (29) 

This linear term goes like an 2 (~ 0.3 for p = 800 MeV) and the integrand is small compared 
to the wall tension E^ and can be neglected. The term which is linear in 8 p is: 


Sp^ 


f 2 u 2 

A 



8 p —2\J 2 p 0 + 


Pu 2 


p 2 p 0 e po 


A 2 

f 2 u 2 


■f3ie po cos (f) 



(30) 


Using again that p Q ~ (/?i//? 2 ) cos (f > 0 we see this simplifies to 

n 2 2 

5 pb) ~ ™ J dz' ^—2 a 7T 2 8pV 2 (cos (f> 0 )^J . (31) 

Since this term is also proportional to arc 2 , 8 p -C p 0 , and the integral of V 2 (cos0 o ) is small, 
we can also neglect this term. The magnitude of the linear term shows how far away we 
are from the exact solution. This information will be used in what follows for the stability 
analysis. 

Now we consider the most important quadratic term. A similar case involving two 
coupled scalar fields was looked at in [^3|] and we will follow the standard procedure presented 
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there closely. If the field configuration is classically stable, the second variation of the 
energy should be a positive differential operator. This means we must solve the following 
Schrodinger-type eigenvalue problem: 


H 


5p 

8(f> 


= U) 


8 p 

8(p 


where H is the operator, 


h = -dh 


A 2 


;U, 


2 f 2 u 2 

and 1 is the 2x2 identity matrix. The potential U is a 2 x 2 matrix with elements: 

u u = P2 (1 + Po) pP ° ~ Pi pp " cos<t> 0 , 

U n = U 2 1 = Pi e Po sill 0 O , 

U22- Pi e p ° cos 4 >„. 


(32) 


(33) 


(34) 

(35) 

(36) 


If the domain wall solution is a stable one, then the operator H is positive-semidehnite. The 
eigenvalue equations can be decoupled by diagonalizing the matrix U. We should note that 
only the potential term has to be diagonalized when looking for negative energy modes, as 
was done in |23|]. The result is 


U± = ^(a± \Ja? - 4 (ab - b 2 - c 2 )) , 
where a, b , and c are defined as 


(37) 


a = P2 (1 +p 0 ) pP % 
b — fi\ e po cos0 o , 
c — (3i e Po sin </> 0 . 


(38) 


The operator H now takes the following form: 


H = 


-dl + 


A 2 


2 f 2 u 




2 „ 2<-7 + 


0 


-dl + U _ 


(39) 


Since we can immediately see that U + > 0 for all z' due to the fact that a > 0, there does 
not exist any negative eigenvalues corresponding to the first equation in this transformed 
basis. Looking at the second eigenvalue equation, we have 


-dl + 


A 2 


2 pu 2 


U- ip„ = uilP 


(40) 


It is a well known theorem of quantum mechanics that there must exist at least one bound- 
state solution to Eq. ((40|). Due to the fact that our domain wall solution should be invariant 
under translations in space z' —> z* + z' Q , there should be a corresponding zero mode in 
the spectrum of Eq. (jd^). In the high density limit, we recover the familiar sine-Gordon 
equation and we know that there is only one bound state in the spectrum of Eq. ([18]). If 
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the exact solution to the equations of motion (|T9|) and (| 20 |) were known, one would expect 
to see a corresponding mode with a vanishing eigenvalue in the spectrum of Eq. (|40|). As 
the density is lowered ( j3\ < /3 2 ) and the saddle point at 0 = tt is still present, one would 
expect the appearance of a mode with a negative eigenvalue corresponding to instability of 
the domain wall. Due to the fact that there still must be a zero mode in the spectrum, the 
zero mode would become the first excited state of Eq. ( |40| ) and the lowest mode would have 
some negative eigenvalue cUq < 0 corresponding to the instability of the domain wall. The 
problem of stability analysis now reduces to determining the eigenvalues corresponding to 
the bound states of Eq. (H). The appearance of an additional bound state in the spectrum 
as the chemical potential is lowered will be the first sign that the system is approaching the 
point of instability. 

Since the solution corresponding to Eq. ( |2l| ) is not the exact solution but does represent 
a path which winds around the barrier at |E| = 0, it is quite possible that the zero mode 
could show up in the spectrum with a small nonzero eigenvalue. It would show up as a 
true zero mode only when the exact solution to the equation of motion is substituted into 
Eqs. (|3|). 

Although the potential [/_ is nontrivial, we will use a variational approach in order to 
determine the upper bounds on Uq and uf. In choosing a trial wave function, we make the 
observation that the potential LL is quite similar to the same potential which arises when 
analyzing the stability of the sine-Gordon soliton, Eq. (|l^). We will pick our normalized 
trial wave function accordingly: 


i>o 



sech(<r^ / ), 


(41) 


with a being the variational parameter. Note that this trial wave function satisfies the 
required boundary conditions 0 a (z' = ±oo) —* 0. For the first excited state, we must pick 
an odd function of z. We will choose the properly normalized function 


/ 3(7 

= y — tanh(crV) sech(oV). 
Applying the variational principle, we must calculate 


-Ei 2 V) = I (-* + ^0^-) !*>}> 


(42) 


(43) 


and minimize this quantity with respect to a to obtain an upper bound E(a) on the energy 
of the nth state. The integral given by Eq. fl43[) must be performed numerically. 

We now have the framework in place in order to test the stability of our domain wall 
configurations. We will assume the perturbative value for the constants / and u as calculated 
in |22[, f 2 = yu 2 /( 87 r 2 ) and u 2 = Setting /1 = 800 MeV and mn s = 100 MeV, we see that 
the ratio of the coefficients is /?i/~ 0.3. In this case, the linear terms (|SU| ) and m are 
small and the solutions given by Eq. ([EL]) are valid approximations to the exact solutions. In 
Fig. H we show the effective Schrodinger potential U e s = X 2 U- /(2/ 2 w 2 ) for p = 800 MeV. 
For the above choice of parameters, the wall tension given by Eq. (p5|) is 


£ (O) (0o, 


2„,2 


Po 


17.48 


ru 

A 


(44) 
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For the trial wave functions given in Eqs. (|4T| ) and (|42|) , the following results for the two 
bound states of Eq. (^Oj) were obtained: 

J 2 u 2 


Eq\<T min 


(ffinin) < -0.016' 


E?\a min ) < +1.163 


A ’ 
pu 2 

A ' 


(45) 


( 2 ) ( 2 ) 

From this, we can see that E\ ^ Eq and both of these quantities are much less the 
wall tension given by Eq. (f|4|). Even though the ground state energy seems to be negative, 
due to the fact that E\ > we can associate this mode with the zero mode. The 
small nonzero eigenvalue is actually an artifact of our approximations. The appearance of 
a negative mode is merely a consequence of the approximate solutions ( |2lD as discussed 
above. This identification can be verified by increasing the chemical potential. When the 
above calculations are repeated as the chemical potential is increased, we see that the energy 
Eq —> 0. This result is expected due to the fact that the |£| field can be integrated out 
as /i increases and the sine-Gordon type theory is recovered. As h is increased, we also 
see that the eigenvalue of the first excited state uq increases towards the maximum value 
of t/ e ff. Eventually, as q is increased further, there is no longer a first excited bound state 
in the spectrum. The calculation of Eq ; and E\ ; was done with the variational functions 
chosen to be different from Eqs. O and fl42|) . The results obtained were the same order 
of magnitude as stated above (|45p. This supports our interpretation of the E^ state as 
the would-be translational mode if exact solutions are known. The magnitudes of the linear 
terms (^) and ([IT]) are approximately the same order of magnitude of Eq , which supports 
our interpretation of Eq 2 ' 1 as the zero mode. 


V. CONCLUSION 

The main goal of this paper was an analysis of the classical stability of U(l) ^ domain walls 
0. Naively one would expect that decreasing /j from n 1 GeV (when the calculations 
are under control i), we inevitably face the situation in which the domain walls become 
unstable objects due to the fast growth of the coefficient /A Q). This naive expectation 
may not necessarily be correct due to the even faster growth of the coefficient @2 (ED, 
which receives contributions from the formation of II molecules as well as perturbative 
contributions. 

What we have actually demonstrated is that the domain wall solution remains classically 
stable down to h — 800 MeV. In order to assess the stability of the domain walls for 
H c ~ H ~ 800 MeV, one must explicitly include the II contribution in the effective potential 
(which is expected to be the dominant contribution at /i ~ /i c |U|). This would hopefully 
lead to /?i //?2 < 1, ensuring the classical stability of the U(1 )a domain walls. Unfortunately, 
due to the lack of information in this region we cannot generalize our results to below 
800 MeV to He- We can argue that the U(1 )a domain walls remain classically stable down 
to h ~ He due to the faster growth of the coefficient compared to (3i as /i is decreased. 
The ratio /?i //?2 ma y be very sensitive to changes in the instanton size distribution or the 
various form factors. This is a difficult problem and unfortunately I do not know how to 
estimate such contributions for a small chemical potential. 
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We should remark here that the stability of the U(1 )a domain wall implies a classical 
stability of the U(1 )a strings ||, which become the edge of the domain walls. It remains to 
be seen whether these (or other) topological defects will have any impact on the physics of 
neutron stars and other compact stellar objects with high core density /i. 
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FIGURES 



|E|/Io 


FIG. 1. The cross section of the effective potential is shown above for /x = 800 MeV. On the 
left half 4> = it is shown and on the right half (f> = 0,2n is shown. 



z 


FIG. 2. This figure shows the effective Schrodinger potential U e s(z') = (A 2 /2 f 2 u 2 )U-(z') as 
a function of the dimensionless coordinate z' where U- is given by Eq. (flTl) . The presence of a 
negative eigenvalue in the spectrum of this potential is indicative of the instability of the classical 
domain wall solution. The effective potential is shown in this figure for /x = 800 MeV. 
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